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ABSTRACT 

We show that most of hot, optically thin accretion disk models which ignore advective 
cooling are not self-consistent. We have found new types of optically thin disk solutions 
where cooling is dominated by radial advection of heat. These new solutions are thermally 
and viscously stable. 

Subject headings: accretion, accretion disks — instabilities 
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1. INTRODUCTION 



In the standard approach to stationary accretion disk models (Shakura & Sunyaev 1973) 
the energy equation takes into account only local, usually radiative, cooling that balances 
the viscous heating. In most cases global heat transport is neglected. This approach is 
justified as long as the disk is geometrically thin because the ratio of the advection term 
representing the global entropy transport to the viscous heating term is of the second order 
in the relative disk thickness H / where H and R are the semi-thickness of the disk and 
the distance from the central object respectively. In some time-dependent calculations of 
thin models (e.g. Taam & Lin 1984; Lasota & Pelat 1991) the advective heat transport is 
taken into account but it represents only a small contribution to the total energy balance. 
In the slim disk approach (Abramowicz, Lasota & Xu 1986; Abramowicz et al 1988; Kato, 
Honma, & Matsumoto 1988; Chen & Taam 1993) terms up to the second order in H/R are 
taken into account so that the advection is part of the scheme together with the (H / R^ - 
order terms in the radial Euler equation. In the high accretion rate regime, when H/R 
is not negligibly small, the cooling of optically thick slim disks is always dominated by 
advection. In the radiation-pressure dominated slim disk models the advective cooling 
rate per unit surface Qadv satisfies the relation (Abramowicz et al 1986) 



where is the viscous heating rate per unit area. 

Recently Narayan & Yi (1994) constructed self-similar, advection dominated accretion 
fiows with Qadv = Q+ and with a sub-Keplerian angular momentum distribution. 

In the last few years, several hot optically thin disk models have been constructed 
following theoretical developments in the description of the microphysical processes which 
may be relevant for the cooling (see e.g., Svensson 1982, 1984). Attention has been focused 
on the electron-positron pair production and annihilation (Kusunose & Takahara 1988, 
1989, 1990; Tritz & Tsuruta 1989; White & Lightman 1989, 1990; Bjornsson & Svensson 
1991a, b, 1992; Kusunose & Mineshige 1992). On the other hand, the relation between 
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optically thick and thin disks has been studied by Wandel & Liang (1991) and Luo & 
Liang (1994) with purely phenomenological "bridging formulas" for the case of effective 
optical depth around unity. The latter authors have concluded that optically thick and 
thin sequences of solutions are always connected in the M (mass accretion rate) versus E 
(surface density) plane. 

In all those models the hot inner regions of the accretion disk are not geometrically thin 
because of the relative inefficiency of the optically thin cooling which makes H/R Cg/vK 
close to unity ( Cg is the adiabatic speed of sound and vk the Keplerian speed). Advection 
is neglected. Furthermore, they are all thermally unstable. This is not just a technical 
difficulty that could be repaired without modifying the basic properties of the proposed 
models. In most cases, the a posteriori calculated advective heat transport is found to 
be the dominant cooling mechanism at high accretion rates. As an example, we show in 
Figure 1 the equilibrium curves on the (M, E)-plane for various models from Kusunose & 
Mineshige (1992) and Luo & Liang (1994); together with the advection dominated lines 
calculated with equation (10) modified to meet their assumptions and units. In Figure la, 
advection is dominant above the heavy dotted line; the heavy dashed line corresponds to 
H/R = 1 (see eq. [7]). Note that equation (10) corresponds to an advection-dominated 
solution independent of the local cooling mechanism (see also Fig. 3). One can see that, 
most of the high temperature solutions are located in the regime where advective cooling is 
dominant. Since it was neglected, the presented models are physically inconsistent. Note 
also that the geometrically thin disk assumption is not satisfied either. This is less serious 
because solutions inconsistent in this sense often possess all qualitative properties of the 
real configurations. In Figure lb, the advection dominant lines are represented by the 
solid, dotted, and dashed ones corresponding to a = 1, 0.1, and 0.01 respectively. It is 
seen that for a < 0.1, the models are not self-consistent. 



2. ADVECTION DOMINATED OPTICALLY THIN SOLUTION 

These difficulties have their source in neglecting the dominant effect of the advective cool- 
ing. Thermally stable optically thin accretion disk models with if/i? < 1 can be con- 
structed when this effect is included. One can get an insight into the properties of optically 
thin advection - dominated accretion disks by assuming that advection and other (local) 
mechanisms cool a Keplerian disk. 

We thus assume the rotational angular velocity is O = Oyt = ^GM/R(R — Rg)^ 
under the pseudo-Newtonian potential $ = —GM/(R — Rg), where M is the mass of the 
central object and Rg is the Schwarzschild radius (Paczyhski & Wiita 1980). Since the 
optical depth is equal to zero the disk is gas-pressure dominated, P = Pg = TZpT/ fx. Here 
p and T are the density and the temperature of the disk respectively, and p is the mean 
molecular weight assumed to be 0.617. Since we are interested in hot disks, opacity will 
be given by electron scattering: k^s = 0.34. We shall assume that radiative cooling is 
provided by optically thin thermal bremmstrahlung with emissivity per unit area, 

Qbrem = 1-24 X 10^^ H evg s"^ cm-^ (2) 

The viscous heating rate per unit area is given by the standard formula 

where M is the mass accretion rate, g = — (2/3)(c/ln 0/c/ln i?), and the factor / = 1 — 
Q{3RG)/^{R){iRG / RY contains the inner boundary conditions. The advection cooling 
rate is taken in a form (see e.g. Chen & Taam 1993): 

Qa,. = = -^T— = ^-e, (4) 

where E = 2H p is the surface density of the disk, S is the specific entropy, and = 
-[(4 - 3/3)/(r3 - 1 )]((/ In T/c/ In i?) + (4 - 3/3)(c/lnE/c/lni?). Here, /3 = P^/P, Tg = 
1 -|- (4 — 3/3)(7 — l)/[/3 + 12(7 — l)(/3 — 1)], and 7 is the ratio of specific heats. For optically 
thin disks /3 = 1. 





Using equations (2-4), the energy equation, = Qadv + Qbrerm can be written as 

Cm^ - 0.361r^/^5^g-^(aE)m + 3.462 x 10"*^a"V ^"^^(aS)^ = 0, (5) 

where q = l-Rcj/R, r = R/Rcj, m = M/Me- Here Me = 4:7rGM/{cKes) is the Eddington 
accretion rate defined as the Eddington luminosity devided by c^. In deriving equation (5) 
we have used the standard equations: 

ly = ^aCsH, (6) 

(H/Rf = {V2/K,,)r-'/'m{aJ:)-'fg-\ = 2r^q\ (7) 

M = -27rRvr^. (8) 

Equation (5) gives the equihbrium relation between rh and E for a given a and for 
any radius. Since it is a quadratic in rh for a given E one expects at most two possible 
values of rh. It is possible to show numerically that a critical E value exists above which 
there is no solution to rh. This corresponds to a single solution mmax^ 

m^ax = 1.7 X 10V-i/2aV/?"'r'- (9) 

One notes that it scales to a^. The relation of mniax(?") for a given a and is shown in 
Figure 2. One can see that for a < 0.1, rhjanx corresponds to a very low mass accretion 
rate. Furthermore, even for rh < rhmax, only a fraction of the disk can be optically thin; 
and the lower the mass accretion rates, the wider the possible optically thin region of the 
accretion disk. 

The solution of equation (5) on the (M, aE)-plane is shown in Figure 3a and b, to- 
gether with the S-shaped sequence corresponding to the optically thick equilibria. Disk 
parameters of M/Mq = 10, r = 5, and a = 0.1 and 0.01 were used. We also assumed 
is unity. For the case of small viscosity, i.e., a < 0.1 ~ 0.2, the optically thin equilibria 
form a sequence disjoined from that formed by the optically thick equilibria. The optically 
thick sequence consists of three branches: (1) gas pressure supported, radiatively cooled 
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one, which is both thermaUy and viscously stable, (2) radiation pressure supported, radia- 
tively cooled one, which is both thermally and viscously unstable, (3) radiation pressure 
supported, advectively cooled one, which is both thermally and viscously stable. The first 
two branches are described by the standard Shakura-Sunyaev solution, the third one is 
given by the slim disk solution. The optically thin sequence consists of two branches: (4) 
gas pressure supported, radiativelly cooled one, which is viscously stable but thermally 
unstable, and (5) gas pressure supported, advectively cooled one, which is both thermally 
and viscously stable. The fact that it is viscously stable should be obvious from its positive 
slope. Thermal stability can be infered from the fact that above that line the heating rate 
is less than the cooling rate but it follows from a criterion properly modified to include 
non-locality of the advection (Abramowicz et al 1994). Similar three stable branches are 
found in the boundary layer numerical calculations of Narayan & Popham (1993) 

The explicit analytical form of this stable equilibrium was not known before. One 
notes from equation (5) that it corresponds to a relation 

m^QMlC^r^l^g'^q-^a^. (10) 

A relation of the form m = f{r)aYl is characteristic of advection dominated a-disk models 
as first noticed by Abramowicz et al (1988; see their eq. [29]). It is also a property of 
advective dominated sub-Keplerian self-similar flows found by Narayan & Yi (1994), which 
is m O.lr-^/^aE for their solutions with small a and 7 = 5/3. In general the form of the 
function / in the m = /aS relation depends on the equation of state and Vr- 

The parameter of can be estimated approximately for the advection dominated 
branch by assuming the correction terms, /, g and q to be unity. Combining equations 
(7), (10), and equation of state, we have T oc and E oc r~"^/^ . Hence, = 1/ (7— 1) — 1/2, 
which is 1 and 5/2 for 7 = 5/3 and 4/3 respecively. 

One should keep in mind that equilibria discussed in this paper were analysed only at a 
given radius. To get correct radial dependences one should solve the full set of conservation 
equations in a similar way as it was done for the optically thick slim disk models. Detailed 
global properties of advection dominated accretion disks will be studied in a future work. 
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FIGURE CAPTIONS 

Figure 1. — (a) Thermal equilibra from Kusunose & Mineshige (1992), a = 0.1. The sohd, 
dashed, long dashed, and dash-dotted lines are for solutions with e~^e~ pairs, without pairs, 
one-temperature, and optically thick disks respectively. Advection is dominant above the 
heavy dotted line. The heavy dashed line correspond to H/R = 1. (b) Thermal equilibra 
from Luo & Liang (1994). The solid triangles, open triangles and stars are for solutions 
with a of 1, 0.1, and 0.01 respectively. The solid, dotted, and dashed lines correspond to 
advection domination limits for the same sequential of a. 

Figure 2. — The variation of rhmax with respect to the radius of the disk. Here = 1 is 
assumed. The solid and dotted lines are for a = 0.1 and 0.01 respectively. Note that for 
a < 0.1, the disk can be optically thin only at very low mass accretion rates. 

Figure 3. — (a) Thermal equilibria for optically thick (the right solid S-shaped line) and 
optically thin (the left solid line) accretion disks. The upper branches represent advection 
dominated solutions. Configuration above the dotted lines r = 1 are optically thin, where 
T is the effective optical depth calculated by assuming that the pressure is dominated either 
by radiation (the upper one) or by gas (the lower one). The parameters assumed here are 
M/Mq = 10, r = 5, a = 0.1, and = 1. (b) The same as (a) except for a = 0.01. 



